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A Possible Origin of Magnetic Fields in Galaxies and 
Clusters: Strong Magnetic fields at z ^ 10? 
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ABSTRACT 

We propose that strong magnetic fields should be generated at shock waves associated 
with formation of galaxies or clusters of galaxies by the Weibel instability, an instability 
in collisionless plasmas. The strength of the magnetic fields generated through this 
mechanism is close to the order of those observed in galaxies or clusters of galaxies 
at present. If the generated fields do not decay rapidly, this indicates that strong 
amplification of magnetic fields after formation of galaxies or clusters of galaxies is 
not required. This mechanism could have worked even at a redshift of ~ 10, and 
therefore the generated magnetic fields may have affected the formation of stars in 
protogalaxies. This model will partially be confirmed by future observations of nearby 
clusters of galaxies. Mechanisms that preserve the magnetic fields for a long time 
without considerable decay are discussed. 
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1 INTRODUCTION 

The question of the origin of galactic magnetic fields is one of 
the most challenging problems in modern astrophysics. It is 
generally assumed that magnetic fields in spiral galaxies are 
amplified an d maintained by a dynamo through rotation of 
the galaxies IIWidrowl2002ll . The dynamo requires seed fields 
to be amplified. However, observations of microgauss fields 
in galaxies at moderate reds hifts strongly constra in the lower 
boundary of the seed fields iAthreva et al.lll9 98h Moreover, 
magnetic fields are also observed in elliptical galaxies and 
galaxy clusters, in which ro tation cannot play a central role 
as the dynamo mechanism ([Clarke. Kronberg. fc Bohringerl 
l200ll: IWidrowlliool IValleell2004D 

The Weibel instability is another m e chanism to gener- 
ate strong magnetic fields fWeibelll95ll lFriedlll959l) . This 
instability is driven in a collisionless plasma, or a tenuous 
ionised gas, by the anisotropy of the particle velocity dis- 
tribution function (PDF) of the plasma. When the PDF 
is anisotropic, currents and then magnetic fields are gen- 
erated in the plasma so that the plasm a particles are de- 
flecte d and the PDF becomes isotropic jMedvedev fc Loebl 
1999). Through the instability, the free energy attributed 
to the PDF anisotropy is transferred to magnetic field en- 
ergy. This instability does not need seed magnetic fields. 
It can be saturated only by nonlinear effects, and thus 



the magnetic fields can be amplified to very high val- 
ues. This instability h as been observe d directly in re- 
cent laser experiments tWei et alj I2002TI . In astrophysical 
plasmas, the instability is expected to develop at shocks 
or at steep temperature gradients, where the PDF is 
anisotropic. Examples of the sites are pulsar winds, shocks 
produced by gamma-ray bursts, jets from active galactic 
nuclei (AGNs), cosmological shocks, and cold fronts (con- 
tact discontinuit ies between cold and hot gas) in clus- 
ters o f galaxies jMedyedev fc Loebl Il999t iKazimura et all 
I l998t iNishikawa et alJ \2Q0A ISchlickeiser fc Shuklal 120031; 
lOkabe fc Hattorill2003h . Although the instability was found 
in 1959, its nonlinear nature had prevented us from under- 
standing its long-term evolution. Recently, however, as com- 
puter power increases, detailed particle simulations of plas- 
mas have been initiated and they have revealed the evolu- 
tion of magnetic fields even after saturat ion of the instability 
JSilva et alJl2003l : lMedvedev et alJl2005lh Based on these re- 
sults, we consider the generation of magnetic fields at galaxy 
and cluster-scale shocks through the Weibel instability at 
the formation of galaxies (both ellipticals and spirals) and 
clusters. We use the cosmological parameters of fio = 0.3, 
A = 0.7, the Hubble constant of Ho = 70kms _1 Mpc -1 , and 
o- 8 = 0.9. 
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2 MODELS 

2.1 Generation of Magnetic Fields at Shocks 

At the vicinity of the shock front of a collisionless shock, 
particles from the upstream are mixed up with those in the 
downstream, and an anisotropy of PDF will be gene r ated m 
the plasma. As discussed by Schlickci ser fc Shuklal J.2003), 
the particles from the upstream region firstly will be affected 
by the Langmuir instability. Since the Langmuir instability 
is a longitudinal electrostatic mode, the velocity component 
parallel to the shock normal will be thermalized, while the 
other components will remain unaffected. Therefore, it is 
natural to assume that the thermal velocity parallel to the 
shock normal is on the order of the relative velocity be- 
tween the upstream and the downstream or on the order 
of the shock velocity Vsh, and those perpendicular to the 
shock normal are on the order of the thermal velocity of the 
upstream plasma. This velocity or temperature anisotropy 
should develop the Weibel instability. 

For shocks in electron-proton plasmas, 
ISchlickeiser fc Shuk la (2003) indicated that the tem- 
perature anisotropy for electrons is too small to drive 
the Weibel instability if M < (m p /m a ) 1/2 = 43, where 
A4 = Vsh/uth.p is the shock Mach number, vth, P is the 
proton thermal velocity in the upstream, and m p and 
m c are the mass of a proton and that of an electron, 
respectively. However, if J\A > 2 (note that M provides 
a measure of the anisotropy in protons), the temperature 
anisotropy for the protons is large enough to drive the 
Weibel instability; this is the case, even if the PDF of the 
electrons is completely isotropic (see Appendix lA"} , 

The magnetic field strength reaches its maximum when 
the Weibel instability saturates, and the saturation level 
would be given as follows. As the Weibel instability devel- 
ops, magn etic fields are generated around n umerous current 
filaments l|Medvedev et al.l2005l:lKatol2005h . The instability 
saturates when the generated magnetic fields eventually in- 
terrupt the current in each filament, or in other words when 
the particle's gyroradii in the excited magnetic fields are 
comp arable to the characteris tic wavelength of the excited 
field jMedvedev fc Loeblll99^ . For electron-positron plas- 
mas, lKatoT|2005}) showed that the typical current at and 
after the saturation is given by the Alfven current (lAlfvenl 
ll939T) . In terms of the Alfven current, the saturated mag- 
netic field strength is given by 

B sat ~ V2I A /(RC), (1) 

where R is the typical radius of a current filament at the 
saturation and Ia = (i>|| iC /c)m e c 3 /e is the Alfven current 
(for nonrelativistic cases). Here, vu e is the electron ther- 
mal velocity in the direction of the higher temperature. For 
strongly anisotr opic c ases (i.e., the particle limit cases dis- 
cussed bv lKato!l2005l) . the radius at the saturation is given 
by R ~ 2c/cj pc , where 0J pe = (47rn e e 2 /m c ) 1 ^ 2 is the electron 
plasma frequency, and the saturated magnetic field strength 
is given by 

B sat -0.5 «|| je (47m c m c ) 1/2 , (2) 

where n c is the electron number density. 

For electron-proton plasmas we consider below, currents 
carried by heavier protons will generate stronger magnetic 
fields at the saturation. It is natural to assume that the sat- 



uration for the proton currents is determined by the Alfven 
current defined for protons, I' A = {v^ p /c)m p c 3 /e, where v\\ p 
is the proton thermal velocity in the direction of the higher 
temperature. In a shock wave, this leads to 

B sat ~ 0.5 V sh (27m p m p ) 1/2 , (3) 

where n p is the proton number density and we used the as- 
sumption of viip ~ Vsh- This expression is valid whe n the 
PDF of the protons has a strong anisotropy or M > 2 iKatol 
2005). Note that the energy density of the saturated mag- 
netic fields attains to a sub-equipartition level with the par- 
ticle kinetic energy. The numbers in the parentheses in equa- 
tions @ and © are different by a factor of two, because in 
electron-positron plasmas, both electron and positron com- 
pose currents at the same time, while in electron-proton 
plasmas, only protons contribute to currents when the pro- 
ton Weibel instability develops. In electron-proton plasmas, 
if electron anisotropy is sufficiently large, electrons first gen- 
erate magnetic fields, which saturate at the early stage of the 
instability, and then protons generate magnetic fields, which 
leads to the second saturation. Su ch an evolution was als o 
observed in numerical simulations iFrederiksen et al"ll2004f) . 

The magnetic fields generated at the shock front will be 
convected downstream at the fluid velocity (in the shock rest 
frame). At least for a short period after the saturation, the 
magnetic field strength decreases as the current filaments 
merge together, because the current in each filament is lim- 
ited to the proton Alfv en current w hile the size of the cur- 
rent filaments increases jKatol2 005). However, the long-term 
evolution is still an open question. iMedvedev et alJ ll2005l) 
indicated that the field should reach an asymptotic value 
beyond which it will not decay, although their model can 
be applied only when currents are not dissipated. Recent 

particle simulations have also shown that the final magnetic 

1/2 

field strength is gi ven by Bf = r? m cr-Bsat, where r/ mcl ~ 0.01 
JSilva et alJl2003t IMedvedev et alJl2005f) . In the following, 
we assume that the magnetic fields of Bf are preserved. The 
very long-ter m evolution of the ma gneti c fields, for which 
the results of lMedvedev et alJ d2005l) and lSilva etafl J2003h 
cannot be directly applied, will be discussed in Section ^] 

We note that unfortunately, at present, most simula- 
tions treat electron-positron shocks and there are virtu- 
ally no simulations that reveal the evolution of collisionless 
electron-proton shocks satisfactorily (in terms of box sizes, 
duration times, and so on), because of the lack of compu- 
tational power. Therefore, the model stated above and used 
in this paper should be confirmed by direct numerical sim- 
ulations in the future. 

2.2 Shock Formation 

According to the standard hierarchical clustering scenario 
of the universe, an initial density fluctuation of dark matter 
in the universe gravitationally grows and collapses; its evo- 
lution can be approxim ated by that of a spherical un iform 
over-dense region llGunn fc Gottlll972l: |Peebleslll980f) . The 
collapsed objects are called 'dark halos' and the gas in these 
objects later forms galaxies or clusters of galaxies. At the 
collapse, the gas is heated by the 'virial shocks' to the virial 
temperature of the dark halo, TUr = GM/ (2r v i r ), where G is 
the gravitational constant, and M and r v i r are the mass and 
the virial radius of the dark halo, respectively. The relation 
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between the virial radius and the virial mass of an object is 
given by 



3M 



47rA c (z)p cr it(.z) 



1/3 



(4) 



where p C iit(z) is the critical density of the universe, and 
A c (z) is the ratio of the average density of the object to the 
critical density at redshift z. The critical density depends on 
redshift because the Hubble constant depends on that, and 
it is given by 



Pcrit(z) 



Pcrit,0^o(l + Z) 3 

Q(z) 



(5) 



where p C rit,o is the critical density at z = 0, and Q(z) is the 
cosmological density parameter given by 



Q(z) = 



n (l + z) 3 
Q {l + z) 3 + X 



(6) 



for the flat universe with non-zero cosmological constant. 
The ratio A c (z) is given by 

A c (z) = 18 7T 2 + 82a: - 39x 2 , (7) 

for the fiat universe l|Brvan fc NormarJll998h . where the pa- 
rameter x is given by x — f2(z) — 1. The virial shocks form 
at r w r v ir and the velocity is v v i r ~ y GM/r v \i. 

In addition, recent cosmological numerical simulations 
have shown that 'large- scale structure (LSS) shocks' form 
even before the collapse dCen fc Ostrikerll99alMiniati et ail 
l2000l: iDave et aljl200ll: iRvu et alJl2003ft . Thev form at the 
turnaround radius (r ta ~ 2r v i r ), the point at which the 
density fluctuation breaks off from the cosmological expan- 
sion. For simplicity, we assume that r ta = 2r v i r , which is 
close to the self-similar infall solution for a particular mass 
shell in the Einste in-de Sitter Universe (rta = 1/0.56 r v i r ; 
lBertschinge3ll985l) . The gas that later forms a galaxy or 
a cluster passes two types of shocks; first, the gas passes 
the outer LSS shock, and then, the inner virial shock. 
The typical velocity of t he LSS shocks is V s h w H(z)r p 
jFurlanetto fc Loeblliooi) . where H(z) is the Hubble con- 
stant at redshift z, and r p is the physical radius that the 
region would have had if it had expanded uniformly with 
the cosmological expansion. The temperature of the post- 
shock gas is T B ~ 3/16(fj,m p /kB)V^ 1 , where fj,m p is the mean 
particle mass, and A?b i s the Boltzmann constant. Note that 
although the model of iFurlanetto fc Loebl |.2004) has been 
compared with numerical simulations at z ~ 0, it has not 
been at high-redshifts; it might have some ambiguity there. 
We do not consider mergers of objects that have already col- 
lapsed as the sites of magnetic field generation because the 
Weibel instability applies only to initially unmagnetised or 
weakly magnetised plasmas; at the merger, collapsed objects 
just bring their magnetic fields to the newly born merged 
object. 

Since the Weibel instability develops in ionised gas 
(plasma), we need to consider the ionisation history of the 
universe. After the entire universe is ionised by stars and/or 
AGNs [z < 8), magnetic fields are first generated at the LSS 
shocks. In this case, we do not consider the subsequent gen- 
eration of magnetic fields at the inner virial shocks, because 
the strength is at most comparable to that of the magnetic 
fields generated at the LSS shocks. On the other hand, when 



the universe is not ionised (z > 8), the Weibel instability 
cannot develop at the outer LSS shocks. However, if the 
LSS shocks heat the gas (mostly hydrogen) to T s > 10 4 K 
and ionise it, the instability can develop at the inner virial 
shocks. 

In this case, the gas ionised at the LSS shock may re- 
combines before it reaches the virial shock. The recombina- 
tion time-scale is given by 

T 



where a is the recombination coefficient, T is the gas tem- 
perature, y is t he ionisation fraction . and uh is the hy- 
drogen density ijShaoiro fc Kan3ll987h . If we assume that 
Tree = Tdyn, where r dy n ~ (l/2)r ta /Vih is the time-scale that 
the gas moves from the LSS shock to the virial shock, the 
ionisation rate when the gas reaches the virial shock is 



Tdyn 



1.22 x 10 s yr 



T 



10 4 K 



mi 



(9) 



for y < 1. We found that y < 1 for z > 9, and the minimum 
value when the generation of magnetic fields is effective (z < 
12, see is y ~ 0.3. When y < 1, we simply replace n p in 
equation with ynu- For temperature, we assumed that 
T — T a in equation Since the magnetic fields do not 
much depend on the temperature (B BS .t oc y 5 oc j l0 35 ) i they 
do not much change even when radiative cooling reduces the 
temperature; at least _B sat does not change by many orders of 
magnitude. It was shown that the ionisation rate just behind 
a shock is y ~ 0.1 , if the shock velocity is relatively sma ll 
(V sh ~ 40 km s-^ IShapiro fc Kandll987tlSusa et alJll998h . 
If the shock velocity is larger, y ~ 1. In our calculations, 
the velocity of the LSS shocks is > 40 km s _1 , when the 
generation of magnetic fields is effective. Thus, y just behind 
the shocks is at least comparable to that obtained through 
the condition of r roc ~ 7"d y n, and the incomplete ionisation 
does not much affect the results shown in the next section 
(see also Section 



3 RESULTS AND DISCUSSION 

FigQ shows the typical mass of objects, M, as a function 
of redshift z; the labels la, 2a, and 3<r indicate the am- 
plitudes of initial density fluctuations in the universe from 
which the objects f orm, on the assumption of the CDM fluc- 
tuations spectrum fearkana fc Loebll200ll) . Fig|5| shows the 
downstream temperature at the virial shock (T v ir) and that 
at the LSS shock (T 3 ) for the objects; T v i r is always larger 
than T a . The ratio T v i r /T s indicates that M > 4 for the 
virial shock. For the LSS shock, M 3> 1, because the gas 
outside the shock is cold. Thus, equation (J^J can be applied 
to both shocks. In Fig[!|] we present the strength of magnetic 
fields (B c ) at a scale of r v i r for the collapsed objects. We as- 
sume that the entire universe is reionised at z = 8. Thus, 
for z > 8, magnetic fields are generated only at the virial 
shocks if T a > 10 4 K. We assume that B c = Bf and plot the 
lines only when T B > 10 4 K. The recombination is effective 
at z > 10 for the 3cr model. On the other hand, for 2 < 8, 
the magnetic fields are generated at the LSS shocks. We 
consider the compression of the fields while the size of the 
gas sphere decreases from r = r ta to r = r v i r , and thus we 
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Figure 1. Typical masses of objects forming at redshift z. The la- 
bels, la, 2a, and 3a (solid, dotted, and dashed lines, respectively) 
indicate the amplitudes of initial density fluctua tions from which 
the o bjects formed; ~ l-3a is the typical value (Barkana & Loch 
1200 ll) . Objects with masses of M < 10 12 M and M > 10 13 M Q 
correspond to galaxies (ellipticals and spirals) and clusters of 
galaxies, respectively. 



assume that B c = 8 2 ^ 3 Bf. Moreover, we plot the lines only 
for Tvii > 2 x 10 K, because below this tempe rature, gas 
infall is suppressed by photoio nisation heating jEfstathioul 
ll992l:lFurlanetto fe Loebll2004l) . 

In Fig|3 the strength of magnetic fields generated by 
protons reaches ~ 10 _8 -10 -7 G and is very close to the 
values o bserved in nearby galaxies and clusters of galaxi es (~ 
1(T 6 G: IClarke et"ai]|200lt Iwidrowll2002l : [Valleell2004) . On 
a galactic scale, the gas sphere may further contract to r <C 
r v i r because of radiative cooling. If magnetic fields are frozen 
in the gas, the strength exceeds ~ 10 -6 G. However, if this 
happens, the magnetic energy exceeds the thermal or kinetic 
energy of the gas. As a result, magnetic reconnection may 
reduce the strength. In fact, the equipartition between the 
magnetic energy density and the thermal or kinetic energ y 
density appears to be held in the Galaxy feeck et al.l fl996). 

The strong magnetic fields shown in Fig. ^indicate that 
strong amplification of magnetic fields, such as dynamo am- 
plification, is not required after formation of the galaxies and 
clusters. This is consistent wi th the observations of galac- 
tic magnetic fields at z > 2 iAthreva et al.l Fl99ST) . Future 
observations of higher-redshift galaxies would discriminate 
between our model and strong dynamo amplification mod- 
els; the latter predict much weaker magnetic fields at higher 
redshifts. Moreover, since the predicted galactic magnetic 
fields are comparable to those at present, they might have 
affected the formation of stars in protogalaxies. Fig[3] also 
shows that our model naturally explains the observational 
fact that the magnetic field strengths of galaxies and galaxy 
clusters fall in a small range (a factor of 10) . Since our model 
predicts that magnetic fields are generated around objects, 
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Figure 2. Temperatures behind the virial shocks (T v ; r ; thick 
lines) and those behind the LSS shocks (T s ; thin lines) for objects 
forming at redshift z. Solid, dotted, and dashed lines correspond 
to la, 2a, and 3a fluctuations, respectively (see FigtTl 




1+Z 

Figure 3. Typical magnetic field strengths of objects forming at 
redshift z. Solid, dotted, and dashed lines correspond to la, 2a, 
and 3a fluctuations, respectively (see FigHt. 



magnetic fields in intergalactic space are not required as the 
seed or origin of galactic magnetic fields. 

Some recent numerical simulations indicated that gas is 
never heated to ~ T v i r for less massive objects because radia- 
tive cooli ng is efficient and the shocks forming at > r v j r are 
unstable (iBirnboim fc Dekelll20o3 iKeres et alJl2005l) . This 
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effect may be important for M < 10 12 M@. If it is correct, 
the generation of magnetic fields is only effective for z < 5 
(the 3cr curve in Fig. 0. However, even for M < 10 12 Mq, 
multiple shocks may form at the inner halo when the col d 
gas reaches there and collides each other dKeres et al .120051) . 
Therefore, magnetic fields may be created there. The details 
of these shocks could be studied by high-resolution numeri- 
cal simulations. 

Although it would be difficult to directly observe the 
generation of magnetic fields through the Weibel instabil- 
ity for distant high-redshift galaxies, it would be easier for 
nearby clusters of galaxies. Since clusters are now growing, 
LSS shocks sho uld be developing outside of the viri al radii 
of the clusters dMiniati et aiJl200fl IRvu et al.ll2003f) . Since 
particles are often accelerated at shocks, the synchrotron 
emission from the accelerated particles could be observed 
with radio telescopes with high se nsitivity at low frequencies 
llKeshet. Waxman. fc Loebl Eooil . The total non-thermal 
luminosity (synchrotron luminosity plus inverse Compton 
scattering with cosmic microwave background [CMB] pho- 
tons) is estimated as 



g fMl 2 
rta/u rn p 2 



(10) 



where e is the acceleration efficiency and / is the 



gas fraction o 


a cluster. If we assume e = 0.03 


Tanimori et al. 


1993 iMuraishi et alJl200fJ) and f = 0.15 


Mohr. Mathiesen. & Evrardlll999t). the maximum value of 



L n t is ~ 10 erg s~ . Since the energy density of magnetic 
fields («b) is smaller than that of the CMB (mcmb), most of 
the non-thermal luminosity (L nt ) is attributed to the inverse 
Compton scattering such as ucMBint/(iiB + ucmb), which 
may have been detected in the hard X-ray band (> 20 keV; 
iFusco-Femiano et alj|2004T) . Thus, the synchrotron radio lu- 
minosity is «Bint/(uB + «cmb) < 10 39 erg s _1 . Some of the 
diffuse radio sources observed in the peri pheral cluster re- 
gions ('radio relics') may be this emission dCovoni fc FerettH 
l2004f) . Since the Weibel instability generates magnetic fields 
on the plane of the shock front, the synchrotron emis- 
sion should be polarised perpendicular to the shock front 
jMedvedev fc Loebl Il999ft. which is actu ally observed for 
some radio relics HlGovoni fc Ferettill2004() . The synchrotron 
emission will tell us the positions of the LSS shocks, if they 
exist. If magnetic fields are generated there, they should 
be observed only downstream of the shock. This may be 
confirmed through Faraday rotation measurements of radio 
sources behind the cluster for both sides of the shock, if 
the coherent length of the fields sufficiently increases (see 
Section 4). 



4 ON THE LONG-TERM EVOLUTION OF 
MAGNETIC FIELDS 

The model proposed here has large uncertainties, especially 
on the very long-term (say Gyr) evolution of magnetic fields. 
The typical time-scale of the Weibel instability until the 
saturation is given by the inverse of the proton plasma fre- 
quency, tjpp = (47rn p e 2 /mp) 1/ ' 2 . However, it is evident that 
this time-scale is much smaller than the cosmological time- 
scale we have considered in this paper. The typical scale- 
length of the saturated fields is given by the proton inertial 



length c/ojpp, and it is also much smaller than a kpc scale. 
Thus, there is a 'missing link' between the generated fields 
and those observed in galaxies and clusters at present. 

If the observed magnetic fields have their origin in the 
Weibel instability, there should be another mechanism that 
takes over the i nstabil i ty. If the saturation and current evo- 
lution model of iKatd ll2005h is applied to proton currents, 
the magnetic field strength should decrease after the satu- 
ration as 



B = B mt {R/RT\ 



(11) 



where R is the radius of a filament, and the filament radius 
at saturation, R, is defined for protons (R ~ 2c/uo pp ). This is 
because the current strength is limited to the proton Alfven 
current I' A while the radius of current filaments increases 
through mergers of the filaments. Since R ~ 10 10 cm for 
n p ~ 10 -5 cm -3 (at r ~ r v i r for a cluster at z ~ 0), the 
magnetic field strength should be reduced by a factor of 
~ 10", if relation JTTJ holds even for an long-term evolution 
to a kpc (~ 10 21 cm) scale field. 

One possibility to overcome this difficulty is that the 
currents are carried in another form at later times. At early 
times of the evolution, the state of the plasma is in a kinetic 
regime in which the characteristic scale of magnetic fields is 
comparable to or smaller than the Larmor radius of parti- 
cles. Probably, at later times, the former would increase via 
current mergers, and would become larger than the latter. 
After this situation is realised, the straight currents gen- 
erated by the Weibel instability, which are limited to the 
Alfven current, might be superseded by a kind of drift cur- 
rents, and the plasma might behave as a MHD fluid. In this 
case, mergers of the cylindrical structures, which were origi- 
nally the current filaments, would occur through the recon- 
nection process as usually considered for MHD fluids, and 
therefore the magnetic field strength would not decrease con- 
siderably unless the diffusion of magnetic fields owing to col- 
lisions between charged particles or collisions with neutrals 
becomes effective. 1 This would allow the magnetic fields to 
survive for a long time. For example, in a MHD fluid, for 
the scale of L ~ 10 10 cm and the temperature of a few 
keV, the dissipation time-scale of magnetic fields through 
the collisions between charged particles is tdiss ~ 10 9 yr 
JSpitzerlll962T) . which is the dynamical time-scale of a clus- 
ter. If current mergers make L larger, tdiss al so becomes 
larger . The numerical simulations performed bv lSilva et al.1 
(2003) showed that magnetic fields do not much decay af- 
ter the saturation. This might reflect the transition to a 
MH D fluid. Moreover ; a M HD inverse cascade mechanism 
(e.g. IVishniac fc Choi 1200 ill could provide another process 
to make larger-scale magnetic fields observed in galaxies and 
clusters (> kpc). Since the ionisation rate, y, is fairly large 
fSection l2.2l l. the ambipolar diffusion of magnetic fields (the 
collisio ns wit h neutrals) can be ignored (e.g. eq. 13-57 in 
EI zeijl97Sft . At any rate, studies about the long-term evo- 
lution of magnetic fields are strongly encouraged. 



We note that 'fast' reconnection (Pctschck 1964) may occur on 
time-scales much shorter than those of the particle collisions. 
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5 CONCLUSIONS 

In this paper, we showed that the Weibel instability can 
generate strong magnetic fields in shocks around galaxies 
and clusters. The strength is comparable to those observed 
in galaxies and clusters at present. The mechanism could 
have worked even at z ~ 10. The results are based on the 
assumption that the magnetic fields generated by the Weibel 
instability are conserved for a long time. The validity of this 
assumption must be confirmed in future studies. 
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It is evident that there is no unstable mode if ^ 0, while 
unstable modes exist if fc^„ > 0. Explicitly, the condition 
for instability is given by 

^a4(a b 2 -l) >0. (A6) 

s 

It should be noted that, in an electron-proton plasma for 
example, even when the electron distribution is completely 
isotropic (a c = 1) at the initial time, unstable modes exist 
if the proton distribution is anisotropic (q p > 1). 

The above results show that even if the Mach number 
of a shock is relatively small (a a w 1), the magnetic fields 
generated by the anisotropy of proton distribution slowly 
grow and re ach the saturation value given by equation 
jKatoll2005h . 



APPENDIX A: DISPERSION RELATION OF 
THE WEIBEL INSTABILITY 

Here, we derive the dispersion relation of the Weibel insta- 
bility in a plasma which consists of some species (or popula- 
tions) of charged particles. In the following, each species is 
denoted by a label 's' (s = e for electron, p for proton, and so 
on), and the mass, charge, and number density of the species 
are denoted by m s , q s and n s , respectively. We assume here 
that each species has a bi-Maxwellian distribution 



if (v) 



(2tt) 3 / 2 c 



■ exp 



2a? 



2a 2 a 2 



(Al) 



where the thermal velocity is given by er s for x or y direc- 
tions, and it is given by a s a s for z direction. 

If the thermal velocity in z direction is larger than the 
other directions 2 , that is, a B ^ 1, the linear dispersion re- 
lation of the Weibel mode, which is relevant to the z com- 
ponent of the current density, is given as 



(cfcf 



i [aaCsZ«s) + al - l] = 0, 



(A2) 



where 



47rn s q| 



(A3) 



and Z(() is the plasma dispersion function defined as 



c 



dz. 



(A4) 



It is shown that the dispersion relation (IA2I can have 
purely positive-imaginary solution of ui, i.e., purely growing 
mode. Since the growth rate becomes zero at the maximum 
wave number of the unstable mode, fc max , we can set u> = 
and QZ(( B ) — at k = fc max in Eq. lA2ll to obtain 



2 Note that in this condition the direction of h igher temper- 
ature is opposite to that of other authors (e.g., IWeibeJ Il959t 
iDavidson et al.lll972l) . Nevertheless, this would be more reason- 
able for anisotropy at shock waves. 
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